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Image Transformation 



Image Transformation 


In many cases it will be useful to transfer the 
image to another domain for better processing or 
applying specific operator or perform further 
analysis. 

From these 

1) Enhancement 

2) Restoration e.g. fixing known distortion 

3) Compression 

4) Analysis 



Frequency Domain: 

The transformation domain of interest is the 
Frequency Domain. 

*Note: transformation === mapping 
T: Spatial domain — > Frequency Domain 

f(x,y) — ► F(v,u) 

The transformations of interest in the context 
of this course are:- 

1) Fourier transform 

2) Discrete cosine transform 


Fourier Transform 


1) Single Dimension and Continuous Function 



f(x) e. 2n i ux dx 


F(u) = T(f (*)) 


/(x) = T - 1 



F(t/) e 2n i Vx du 


2) Two Dimensions Continuous function: 

F(u,v) = T(f(x,y)) - f(x, y)e ~ 27Tj( ’ ,x+vy) dxdy 


f{x,y) = :F- 1 (F(U, V)) = fC rn F(U,V)e 27Tj(ux+VY) dvdi> 


Fourier Transform 


3) Single Discrete: 


F(u) 


Given that f{x) = {/( 0),/(l), - 1)) 

Or x = 0,1|2, ... ... ,N — 1 

4) Two Dimension Discrete: 

f(x,y) 
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y = ,m-i 
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Fourier Transform 


4) Two Dimension Discrete: 


F(U,V)= T(f(x,y)) = 
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f(x,y ) is a real domain matrix , F(u f v) is a complex domain 
Signal magnitude: (any element in the matrix have its own magnitude) 


\\F{U,V)\\ = y/R^U.V) + I 2 (U,V) 

, . , I(U,V) 

the phase of an element is <p(F(Lf,F)) = tan 1 ^ 


||F(t/, V)\\ 2 = Signal Power = R 2 (U,V) + I 2 (U,V) 


— d2. 


Discrete Cosine Transform 


Transform of the discrete cosine is from 
Real Domain — > Real Domain 


Single Dimension: 


JV-l 


c(u) = a (u) ^ f(x) 
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c(u) -* f requency Domain 
Two Dimensional Discrete Cosine Transform: 
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Image Enhancement in Frequency 

Domain 


Image Enhancement in Frequency Domain 


/ O, y) 


E nil an c e me n t 


Orignal Image 


► fix, y) 


Enhanced 

Image 


^(F(x,y)) 


F(u t v) 


E nliancemeiit 


all in Frequency Domain 


1) We convert the image to frequency domain 

2) We enhance the image in frequency domain 

3) Then convert it to/(x, y) 


F r (u, i?) 


Discrete Fourier Transform 


• Fourier analysis is fundamentally a method for expressing a function as a 

sum of periodic components, and for recovering the function from those 
components, (DFT) for discrete input function . 

• The DFT has become a mainstay of numerical computing in part because of a 
very fast algorithm for computing it, called the Fast Fourier Transform (FFT). 

• DFT has a great number of applications in digital signal processing, e.g., for 
filtering. 

• the discretized input to the transform is referred to as a signal in time domain. 
The output is called a *spectrum* or transform* and exists in the *frequency 
domain*. 

• Some terms 

- The least Frequency element (zero frequency component will be at top 
left corner) 

- The Highest Frequency element 


Image Processing in the Fourier 

Domain 


the Fourier Transform tells you what is happening in the image in terms 
the frequencies of those sinusoids: 

For example: 


A picture of a plain wall : the 

values of the pixels change very little 
as you go from left to right or from 
top to bottom. In the frequency 
domain that means that your image 
contains low frequencies, but no high 
frequencies. 





Image Processing in the Fourier 

Domain 

A picture of a picket Fence : 

the values of the pixels change all 
the time as you go from left to 
right. So in Fourier domain you 
have high frequencies in the X 
direction, but not in the Y direction. 


input Image Magnitude Spectrum 



A picture of a checkerboard : 

then the pixel values change a lot 
in both directions. Thus the Fourier 
transform of the image will have 
high frequencies in both X and Y. 









Image Processing in the Fourier 

Domain 

Magnitude of the FT 




Convolution i 






i Multiplication in 
r Domain 
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Low-pass Filtering 


Original image 


FFT of original image 



Low-pass filter 




Let the low frequencies pass and eliminating the high frequencies. 


Low-pass image 



FFT of low-pass image 



Generates image with overall 
shading, but not much detail 






High-pass Filtering 


Original image 


FFT of original image 


u'itmri: 


Lets through the high frequencies (the detail), but eliminates the low 
frequencies (the overall shape). It acts like an edge enhancer. 

High-pass image FFT of high-pass image 




Discrete Fourier Transform in 

Numpy 

np.fft.fft2() : an FFT package that provides us the 
frequency transform which will be a complex array. 

— first argument is the input image in gray scale . 

— Second argument is the size of output array.(Optional) 

• If size > size of input image -> input image is padded with zeros before calculation 
of FFT. 

• If size < size input image, input image will be cropped. 

• If no arguments passed, Output array size will be same as input. 

np.fft.fftshift(A) :shifts transforms and their frequencies to 
put the zero-frequency components in the middle, and 
np.fft.ifftshift(A) undoes that shift. 

np.fft.ifftshiftQ & np.fft.ifft2( ) 



Discrete Fourier Transform in 

OpenCV 

cv2.dft() and cv2.idft() : 

— It returns the same result as previous, but with two channels. First channel will 
have the real part of the result and second channel will have the imaginary 
part of the result. 

- The input image should be converted to np.float32 first. 


dft = cv2. dft ( np. float 32 (img)jf lags = cv2.,DFT_COMPLEX_OUTPUT) 
d-Pt_shi-Ft = np.fft.fftshift(dft) 

magnitude_spectrum = 26*np. log (cv2. magnitude (dft_shi-ft [ \ 3 : ^8] j dft_shift[ : ^ : ,1] ) ) 



»> help (np . f f t . f f tshif t) 

Help on function fft3hift in module numpy. f ft .helper: 
f f t3hif t (x, SKe3=None) 

Shift the zero- frequency component to the center of the spectrum. 

This function swaps half-spaces for all axes listed (defaults to all) . 
Note that "'ylO]'' is the Nyquist component only if k 1 len (x) 1 k is even. 

Faram.eters 


x : array_like 
Input array. 

axes : int or shape tuple , optional 

Axes over which to shift. Default is None, which shifts all axes. 

Returns 


y : ndarray 
Examples 


»> freqs = 

np . f f t . f f tf req 

(10, 0.1) 





»> freq3 
array (I 0 . , 

1., 2., 3., 

4., -5., 

-4., 

-3., 

-2., 

- 1 . 1 ) 

»> np.fft. 
array ( [-5 . , 

fftshift (freqs) 
-4., -3., -2., 

1 

M 

Q 

i.. 

2 ., 

3., 

4 . 1 ) 


Shift the zero- frequency component only along the second axis: 

>» freqs = np . f f t . f f tf req ( 9 , d=i . /9 ). reshape ( 3 , 3) 

»> freqs 

array ( [ [ 0., 1. , 2.], 

I 3., 4., -4.], 

1-3., -2., -1.11) 

»> np . f ft . fftshift (freqs, axes=(l,)} 


array ( I [ 2 . , 

o.. 

1-1, 

1-4., 

3., 

4.1, 

[-1., 

-3., 

-2.11) 


